Heavy Quarkonium Production at LHC through W Boson Decays 
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The production of tiie iieavy (cc)-quarkonium, (c&)-quarkonium and (66)-quarkonium states 
[{QQ') quarkonium for short], via the semi-inclusive decays, has been systematically studied 
within the framework of the non-relativistic QCD. In addition to the two color-singlet S-wave 
states, we also discuss the production of the four color-singlet P-wave states |(QQ')(^'Pi)i) and 
\{QQ')C Pj)i) [with J = (0, 1,2)] together with the two color-octet components \ (QQ')C So)s) and 
\{QQ'){^ Si)s) ■ Improved trace technology is adopted to derive the simplified analytic expressions 
at the amplitude level, which shall be useful for dealing with the following cascade decay channels. 
At the LHC with the luminosity £. oc 10 cm~ s~ and the center-of-mass energy ^/S = 14 TeV, 
sizable heavy-quarkonium events can be produced through the VK"*" boson decays, i.e. 2.57 x 10*" rjc, 
2.65 X 10® J/* and 2.40 x lO*^ P-wave charmonium events per year can be obtained; and 1.01 x 10 
Pc, 9.11 X lO'* P* and 3.16 x lO" P-wave (c&)-quarkonium events per year can be obtained. Main 
theoretical uncertainties have also been discussed. By adding the uncertainties caused by the quark 
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I. INTRODUCTION 

The study of the W boson is helpful for understand- 
ing the electrowcak interactions and for searching new 
physics beyond the standard model. With the LHC lu- 
minosity rising up to £ cx 10^'^cm~^s~^ and running at 
the center-of-mass energy y/S = 14 TeV, large amount of 
W bosons about 10^" events per year will be produced. 
This makes the LHC a much better W boson factory than 
the TEVATRON and more W boson rare decays 

can be adopted for precise studies. In Refs.[l,|^, the au- 
thors have discussed a class of W boson semi-inclusive 
decays to the lowest S'-wave heavy-quarkonium states 
{QQ') with Q and Q' stands for the c-quark or the b- 
quark respectively. Their results show that large number 
of heavy-quarkonium events through the W boson decays 
can be found at LHC, so these channels shall be helpful 
for studying heavy-quarkonium properties. 

Intuitively, the heavy-quarkonium production process 
could be understood in terms of two distinct steps: the 
production of the QQ' pair and the subsequent evolution 
of the QQ' pair into the quarkonium. Different treat- 
ment of the evolution leads to different theoretical mod- 
els, among which the non-relativistic QCD (NRQCD) 
is widely adopted. In the framework of NRQCD, a 
doubly heavy meson is considered as an expansion of 
various Fock states. And in addition to the two color- 
singlet S'-wave states \{QQ')i^ So)i) and \{QQ'){^ Si)i) , 
the naive NRQCD scaling rule shows that the four color- 
singlet P-wave states \iQQ'){^ Pi)i) and {QQ')i^Pj)i) 
[with J = (0, 1, 2)] together with the two color-octet com- 
ponents \{QQ')CSo)s) and \{QQ')i^Si)s) shall also give 
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sizable contributions to the production. Here the thick- 
ened subscripts of {QQ') stand for the color indices, 1 for 
color singlet and 8 for color octet; the relevant angular 
momentum quantum numbers are shown in the parenthe- 
ses accordingly. These higher excited (QQ'j-quarkonium 
states may directly or indirectly decay to their ground 
state via the electromagnetic or hadronic interactions 
with high probability. It is interesting to study higher 
Fock states' contributions to make a sound estimation 
on the heavy-quarkonium production, and hence to be a 
more useful reference for experimental studies. 

Moreover, the heavy-quarkonium production itself is 
very useful for testing perturbative QCD ilrQ. For ex- 
ample, since its discovery by the CDF collaboration [lo| . 
the Be meson being the unique 'doubly heavy-flavored' 
meson in the standard model has aroused people's great 
interest. The 'direct' hadronic production of the Be me- 
son has been studied systematically in Refs. (Tll - [l^ . As 
a compensation, it would be helpful to study its 'indi- 
rect' production mechanisms. Because sizable top-quark 
and W boson events shall be produced at the LHC, the 
production of Be through their decay shall be helpful for 
determining the i3c-meson properties, since too many di- 
rectly produced Be events shall be cut off by the trigging 
condition at the LHC 0-Ill • A systematical study on the 
Sc-meson production through the top-quark and the Z° 
boson decay can be found in the literature fl5l - l20j . In 
the present paper, we shall make a systematic study on 
the Be meson production through the W boson decays. 

To deal with the heavy-quarkonium production 
through the W boson semi-inclusive decays, one needs 
to derive the squared amplitude, which is usually done 
by the conventional trace technique. The analytical ex- 
pression for the squared amplitude of the S-wave case can 
be found in Ref . [J] , however it is hard to write down the 
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squared amplitudes for the P-wave cases, which is much 
too complex and lengthy. One important way to solve 
this is to deal with the process directly at the amplitude 
level. For this purpose, the helicity amplitude approach 
and the improved trace amplitude approach have been 
suggested in the literature. As for the helicity amplitude 
approach [2lt , all the amplitudes can be expressed by the 
complex valued helicity amplitudes that can be numer- 
ically calculated, an explicit example of which to deal 
with the case of massive spinors can be found in Ref . [3] . 
While for the improved trace amplitude approach sug- 
gested and developed by Refs. pi5l - [l9| . the hard-scattering 
amplitude can also be expressed by the dot-products of 
the particle momenta as that of the squared amplitude, 
it is, however, done at the amplitude level and is much 
more simpler. In the present paper, we shall adopt the 
improved trace amplitude approach to derive analytical 
expressions for all the mentioned Fock states, and to be 
a useful reference, we simplify its form as compactly as 
possible when fully applying the symmetries and rela- 
tions among them. 

The paper is organized as follows. In Sec. II, we show 
our calculation techniques for the mentioned semi- 
inclusive decays to the heavy-quarkonium. In Sec. Ill, we 
present the numerical results and discuss on the prop- 
erties of the heavy-quarkonium production through W'^ 
decays. The final section is reserved for a summary. To 
make the paper more compact, we present the detailed 
formulas for dealing with the process under the improved 
trace amplitude approach in the Appendix. 

II. CALCULATION TECHNOLOGY 

We shall deal with some typical W boson semi- 
inclusive processes for the heavy-quarkonium production, 
i.e. W+{k) {QQ')[n]{q3) + qiq2) + q'iqi), where q and 
q' stand for the Cabibbo-Kobayashi-Maskawa (CKM) fa- 
vored quark and anti-quark accordingly, and k and qi 
are momenta of the corresponding particles. According 
to the NRQCD factorization formula [l^, its total decay 
width dT can be factorized as 

dr = Y, dr{W+ ^ {QQ')[n] + qq') (0« (n)) , (1) 

n 

where {O^ (n)) describes the hadronization of a QQ' pair 
into the observable quark state H and is proportional 
to the transition probability of the perturbative state 
{QQ')[n] into the bound state. As for the color-singlet 
components, their matrix elements can be directly re- 
lated to the wave functions at the origin for the S'-wave 
states or the first derivative of the wave functions at 
the origin for the P-wavc states [d], which can be com- 
puted via the potential models [23l428l | and/or potential 
NRQCD m, [S^ and/or lattice QCD [HI respectively. 
As for the color-octet components, their matrix elements 
are to be determined experimentally, which are smaller 
than the color-singlet matrix elements by a certain v"^ 



order. Here v is the typical velocity of the heavy quark 
or anti-quark in the quarkonium rest frame, v'^ ~ 0.3 for 
J/^* and ~ 0.1 for T. More specifically, based on the 
velocity scaling rule and under the vacuum-saturation 
approximation, we have H, 0, [H, IHj 

{{QQ%['So]\OsCSo)\{QQ'M'So]} 
~ As{v)^{{QQ')i['So]\OiCSo)\{QQ')i['So]) (2) 

and 

{{QQ')s['Si]\OsCS,)\{QQ%['Si]) 
^ As{v)^{{QQ')i['S^]\0,eS,)\{QQ'U['S^]) , (3) 

where A5(w) is of order v^. 

The short-distance decay width 

dr{W+ ^ {QQ')[n] + qq') = ^^jMpd^g, (4) 

where ^ means that we need to average over the spin 
states of the initial particles and to sum over the color 
and spin of all the final particles. In the rest frame, 
the three-particle phase space can be written as 

The 1 — >■ 3 phase space with massive quark/antiqark in 
the final state can be found in Refs.flg. [l8| . To shorten 
the paper, we shall not present it here. With the help 
of the formulas listed in Refs.[l6l. [T8| . one can not only 
derive the whole decay width but also obtain the corre- 
sponding differential decay widths that are helpful for ex- 
perimental studies, such as dT/dsi, dT/ds2, dT/dcos9i3 
and dr/dcos923, where Si = (gi -I- 93)^, S2 = (gi -I- 92)^, 
6'i3 is the angle between qi and ga, and 823 is the angle 
between g2 and ga in the rest frame, respectively. 

And then our task left to deal with is the hard- 
scattering amplitude for the specified processes 

W+ {cc)[n] + cs, W+ {hh)[n] + cb, 
W+ -> {ch)[n] + bs, W+ -> {cb)[n] + cc. 

Their amplitudes can be generally expressed as 

m 

iM = Cusi[q2)'^AnVs'j[qi), (6) 

where m stands for the number of the Feynman dia- 
grams, s and s' arc spin indices, i and j are color indices 
for the outgoing quark and antiquark. The overall fac- 
tor C = Cs or Co stands for the specified quarkonium in 
the color-singlet and the color-octet states respectively. 
= and Co - ^^^^(V2r"T^r«),„ 

where y/2T^ stands for the color factor of the color- 
octet quarkonium state. Vckm stands for the CKM ma- 
trix element, Vckm = for — > (cc)[n] + cs and 
W+ {d))[n] + bs- Vckm = Kb for W+ {cb)[n] + cc 
and W+ (bb) + cb. 
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FIG. 1. Feynman diagrams for the precess W^{k) — > 
{cQ')(qs) + Q'(g2)s(gi), where {cQ') stands for the quarko- 
nium Fock states \{cQ')i[\So]) , |(cQ')i P^i]) , |(cQ')i['A]), 
\{cQ')s['So]g), l(cQ')iPf'j]> and \{cQ')s\'' Si]g) respectively. 



A. An for VK+ (cc)[n] + cs and -> {cb)[n] + bs 

For convenience, we present these two processes as 
W+{k) (cQ')N('?3) + Q'{q2)s{qi), where Q' stands 
for c or b quark accordingly. The Feynman diagrams of 
the process arc presented in Fig.([T|), where the intermedi- 
ate gluon should be hard enough to produce a cc pair or 
bb pair, so the amplitude is perturbative QCD calculable. 

For {cQ') quarkonium in S'-wave states. An can be 
written as 



A2 



" (932 + q2Y 



(qi + 92 + 932)^ - m 



J 9=0 

la 



J 9=0 



And for the P-wave states, An can be written as 



,s=o,L=i m \ " 
d 



A. 



S=0.L = 1 \ 



7a 



7a 



('732+92)^ (g2 + 93)^-nT-c 



(932 +92)^ 



((71 + 92 + 932)^ - "l2 



J 9=0 

7a 



9=0 



and 



.S=1,L=1 _ \ " 



■ ^2 + ^3 

'a / , NO Ja 



A 



S=\,L=\ _ J 



dqfj 



7a 



(932 + 92)^ (92 + 93)^ - ml 

(932 + 92) 



9=0 



(91 + 92 + 932)^ - m; 



r7a 



9=0 



(7) 
(8) 

(9) 
(10) 

(11) 
(12) 



Here q stands for the relative momentum between the 
two constituent quarks in {cQ') quarkonium. q^i and 532 
are the momenta of the two constituent quarks, i.e. 



931 



M 



93 + 9, 932 



M 



93 



(13) 



where M ~ to^ + "t-q' are adopted to ensure the gauge 
invariance of the hard-scattering amplitude. e{k) is the 
polarization vector of W^. Ssiqa) and e;(93) are the po- 
larization vectors relating to the spin and the orbit an- 
gular momentum of (cQ') quarkonium, and £^1^(93) is 
the polarization tensor for the spin-triplet P-wavc states 
with J = 0, 1 and 2 respectively. The projectors 11°^ (g) 
and (g) are for spin-singlet and spin-triplet quarko- 
nium states respectively, and their covariant form can be 
conveniently written as 



n°(9) 



M 



"32 



"iQ')75(^3^ +n^c) (14) 



and 



AnicmQi 



"32 



-"iQ')7t'(^3i +"T-c)- (15) 



After substituting these projectors into the amplitudes, 
the amplitudes then can be squared, summed over the 
freedoms in the final state and averaged over the ones in 
the initial state. Selection of the appropriate total angu- 
lar momentum quantum number is done by performing 
the proper polarization sum, which for a spin-triplet S 
state or a spin-singlet P state is given by 



E 



where 11 



a/3 



-5a/3 



and Jz 



(16) 



Sz or 



respec- 



tively. And for the case of 3pj states, the sum over po- 
larization is given by (22l . [33j 



^(0)^(0)* 
=-a^^a',3' 
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FIG. 2. Feynman diagrams for the process VK^(fc) — > 
{Qb){q3) + c(f?2)Q(gi), where (cQ') stands for the quarko- 
nium Fock states \{Qb)i[' So]) , \{Qb)^['Pi]), 
\iQb)s['So]g), \{Qb)ifPj]) and \{Qb)s\:'Si]g) respectively. 



A2 = 



And for the P-wave states, An c£ 

aS=0,L=1 _ u 

"^1 — 



.5=0, L=l _ n 

and 

^S=l,L=l _ J 
.S=1,L=1 _ 7 

The quark momenta q^i and (732, and the projectors 
n"^ (q) and 11^^ (q) have the same meaning as those hsted 
in the last subsection, and one needs to change the two 
constituent quarks there (c and Q') to be (Q and b). 



C. Analytical expressions for the mentioned 
amplitudes 

As mentioned in the Introduction, we adopt the im- 
proved trace technology to deal with the hard-scattering 
amplitude ([5]) . Detailed processes of the approach can be 
found in Refs. [l6l - [l9| . here for self-consistency, we shall 
present its main idea and our main results. 

For the present considered W boson semi-inclusive de- 



E(l) (1)* _ ^aa'^PP' — ^aP'^a'P 

and 

E(2) (2)* _ n„a'n^^/ -f IlQ^'IIa/^ ^aP^a' P' 
^aP^a'P' ~ 2 3 



B. An for W+ (cb)[n\ + cc and W+ ibb)[n] + cb 

For convenience, we present these two processes as 
W^{k) ^ {Qh)[n]{q3,) c{q2)Q{qi), where Q stands for 
the c or b quark accordingly. The Feynman diagrams 
of the process are presented in Fig.®. For the {Qb) 
quarkonium in S'-wave states. An reads 



(17) 



(18) 



(19) 



(20) 



(21) 
(22) 

I 

cay process, W+{k) iQQ')[n]{q3) + q{q2) + q' (qi), there 
are one quark and one antiquark in the final state. Un- 
der the approach, we first arrange the whole amplitude 
into four orthogonal sub-amplitudes M±s±s' according 
to the spins of the outgoing quark with spin-index s and 
antiquark with spin- index s', then transform these sub- 
amplitudes into the trace form by properly dealing with 
the massive spinors with the help of an arbitrary light- 
like momentum ko and an arbitrary space-like momen- 
tum ki, ki = —1 and ko ■ ki — 0. And then we do the 
trace of the Dirac 7 matrix strings at the amplitude level. 



The final results are independent of kg and ki, which provides 
another way to check the rightness of the derived results. And 



5^ "^f-^i 



^3 



(qi +93)^ -ml^°'{q3i +gi)2^" 



la 



r^(fc)(l-7' 



[qi + 92 + 931)^ - 



('731 + qi) 



2 7a 



J 9=0 



m be written as 
d 



(93) 



dq^ 



{qi + qzY - ml (ggi +qi) 



q=0 



(.3);^ 



la 



'1 ^2 T31 ^.^(fc)(l„^5) "93 



[qi + q2 + q?,iY - ml 



(qai + qi)' 



:1a 



J 9=0 



.(93) 
.(93) 



d 

dqf, 
d 
dqt. 



{qi + q^Y-ml [qzi+qiY 



9=0 



la 



^1 + ^2 + ^31 + 

{qi + (72 + 931)^ - m\ 



r^(fc)(l-7') 



(931 + 91)' 



r7a 



9=0 
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which shall result in explicit series over some independent 
Lorentz-structures. 

After summing up the spin states of the outgoing 
quark/antiquark, the squared amplitude can be divided 
into four orthogonal parts, 

|M|2 = |Mi|2 + IM2P + |A/3|2 + IM4P, (23) 

where the four amplitudes Mi can be written as 



and 



M^' [n) = 4 (") = 3, 4) 



(32) 



where m stands for the number of basic Lorentz 
structures Bjin), whose values depend on the heavy- 
quarkonium state [n]. The independent lorentz struc- 
tures Bj [n) for all the Fock states and the explicit expres- 



Ah = 
A/2 = 

M3 = 
and 

A'U = 



Mss' + M-s-s' 
^/2 



NTt 



(^1 -mgOl^o(^2 +"^9)^ 



sions for the Lorentz-invariant coefficients A^-''^{n) and 



A^^ (n) are put in the Appendix. 



Mss' - M^s-s' 



NTt 



(^1 -"i,,)75|^o(^2 +mg)A 



Ms 



M-s 



NTt 



V2 

(^^ - mg')|^o|^i(^2 + "^9)^ 



^/2 



Ms 



NTt 



V2 



^/2 



(24) 
2 

Here, A = ^ ^„ and = 1/^4(^0 • gi)(fco • 92) is 

n=l 

the normalization constant. As a viable choice to sim- 
plify the amplitude, we take fco = (72 ~ oiqi with a = 

^^-^2 ^ and fcj^ = tNae'^''P''qi^kpq2a ( A^o en- 
sures fci • fci 



III. NUMERICAL RESULTS 

We adopt the following values to do the numerical cal- 
culation Hill: mw = 80.399GeV, = 2.085GeV, 
rric = 1.35GeV, = 4.90GeV, = O.lOSGeV, \Vcs\ = 
1.023± 0.036 and \Vcb\ = 0.0406± 0.0013. Leading-order 
as running is adopted and we set the rcnormalization 
scale to be 2mc for charmonium and {cb) quarkonium, 
■and 2mf, for bottomonium accordingly, which lead to 
Q!s(2toc) = 0.26 and as{2mi,) = 0.18. Non-perturbative 
matrix elements can be related to the wave function at 
the origin 5's(0) = y/l/ATTRg{0) and the first derivative 
of the wave function at the origin ^'p(O) = ^/3/A^TR'p{0), 
where we adopt 



-1), which leads to 



i^k ■ q-z-qi- q2$ - 



|i?s(cc)(0)p 0.810 GeV^ , \R'p{cc){0)\ 
|i?s(c6)(0)p = 1.642 GeV^ , |i?'p(c5)(0)p 
|i?s(6&)(0)|2 = 6.477 GeV^ , |i?^(66)(0)p 



2 0.075 GeV^ 
0.201 GeV^ 



1.417 GeV'^ 



Then the resultant Ali are. 

Ml = ii X Tr[{<^^ ~ m,0(^2 + 
A'/2 = i2 X Tr[{(^^ - m,075(^2 + ™'/)^] 
Af3 - My - N^[mq{qi ■ k) + m,,{q2 ■ k)]M2 (27) 
Mi = Mi, + Na[mg{qi ■ k) - m,-(g2 • k)]Mi (28) 



(25) 
(26) 



where Li^2 = l/(2^gi • (?2 T ' 



') and 



M3, = ^Tr 

4i2 

A^o 

A/. . -^T. 



(^1 -™g')75l^(^2+'^9)^ 
2 + ?71q)^ 



(29) 
(30) 



As a cross-check, in addition to the improved trace 
technology, we also adopt the traditional trace technol- 
ogy for dealing with the mentioned processes. Numeri- 
cally, we obtain a nice agreement between these two ap- 
proaches for all the above mentioned decay channels and 
heavy-quarkonium states. Moreover, it is found that our 
numerical results for the color-singlet 5-wave cases agree 
with those of Ref.[3l under the same input values. 



A. Basic results 

As a reference, we calculate the decay widths for the 
basic processes c + s and — > c + 6. Their 

decay width can be written as 



Furthermore, the amplitudes Af-o can be expanded 
over some basic Lorentz structures: 



r 



Gf\VckmY 
71 



3^(m2 + |p|2)(m2 + |p|2) + |p|2 



M,{n) = Y,A]{n)B,{n){i = l-A) 



(31) where mi = rric, TO2 = 'm's or mi,, and 



Vim^ - {mi - m2Y){m^ - (mi -t- 7712)^) 
2mvF 



one can choose them to be those that can maximumly simply the 
analytical expressions for the amplitude. 



VcKM = Vcs for W+ cs and Vckm = Vcb for W+ 
cb respectively. Then, we obtain Fyiz+^c+s — 713.6 MeV 
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Vv — ^ (cc)[n] + cs 


r(KeV) 


^ W+->(cc)[n] 


W+ Vc 


174.8 


2.45 X 10"'' 


W+ j/^p 


180.6 


2.53 X 10"'' 


H/+^l(cc)i[iPi]) 


37.9 


5.31 X lO^'^ 


W+^|(cc)ipPo]) 


42.5 


5.95 X 10^5 


H/+^|(cc)ipPi]) 


45.1 


6.32 X lO^'^ 


H/+->!(cc)ipP2]) 


39.9 


5.59 X 10-^ 


W+ -> |(cc)8[i5o]) 


21.9«'' 


3.071)-' X 10^5 


W+ ^ |(cc)spSi]) 


22. 6u'' 


3.17?;'' X 10~5 



TABLE I. Decay widths and branching fractions for the char- 
monium production through — (cc)[n] + cs. 



W+ {cb)[n]+bs 


r(KeV) 


' w+^(cE)H 


W+ 


6.32 


8.86 X 10^^ 


W+ ^ B* 


5.38 


7.54 X lO-'' 


|(c&)i[iPi]) 


0.300 


4.20 X lO"'^ 




0.851 


1.19 X 10"'' 


|(c5)ipPi]) 


0.583 


8.17 X lO"'^ 


M/^ |(c6)ipP2]) 


0.0326 


4.57 X 10"* 


W ^\{cb)s['S»]) 


0.789u'' 


l.lli)-* X lO"'' 


|(cb)8pS'i]) 


0.672U'' 


0.94?;-* X lO"'' 



TABLE IL Decay widths and branching fractions for the (cb)- 
quarkonium production through — > (cb) [n] + bs. 



W+ {cb)[n]+cc 


r(KeV) 


1 w^CcSjln] 


W+ Be 


0.546 


4.88 X 10"" 


W+ B* 


0.810 


7.25 X 10"" 


|(c6)i[*Pi]) 


0.170 


1.52 X 10"" 


|(c&)ipPo]) 


0.037 


3.31 X 10"5 


|(c6)ipPi]) 


0.079 


7.02 X 10"5 


|(c&)ipP2]) 


0.094 


8.36 X 10-5 


|(c6)8[*So]) 


0.0687;" 


6.11i;" X 10-5 


W ^ \{cb)s[' S{\) 


ClOlt;-* 


9.03t;" X 10-5 



TABLE III. Decay widths and branching fractions for the 
(c6)-quarkonium production through — >■ (c&)[n] + cc. 



W+ (66) [n] +CC 


r(eV) 


iV^(66)[„] 


W+ ^ ?76 


17.6 


1.57 X 10"5 


H/+ T 


18.6 


1.66 X 10"5 


W->!(66)i["Pi]) 


0.469 


4.19 X 10""^ 


|(66)ipPol) 


0.771 


6.90 X lO"'^ 


W^^|(66)ipPi]) 


0.821 


7.34 X 10""^ 


|(66)ipP2l) 


0.284 


2.54 X 10-'^ 


1(66)8 ["S-o]) 


2.20t;" 


1.97«" X IQ-'' 


W^-^ 1(66)8 pS'i]) 


2.33?;" 


2.08w" X 10-'' 



TABLE IV. Decay widths and branching fractions for the 
bottomonium production through — >■ (66) [n] + cc. 



and T^+^^j^i = 1.118 McV, the sum of which is about 
34% for the total width Tw+ . 

Total decay widths and their branching fractions for 
the typical channels of {QQ')[n] + q + q' are 

listed in TABs. piIIIIIIIIV|) . It is found that the squared 
amplitude for the color-octet decay width is suppressed 
by eight times to that of the color-singlet case. As a 
combined effect of such color suppression and the rela- 
tive velocity suppression (w^-suppression), the color-octet 
channels are quite small in comparison to their corre- 
sponding color-singlet production channels as shown by 
TABs. piIIIIIIIIV|) explicitly. So in the following discus- 
sion, if not specially stated, we shall not include the color- 
octet states' contributions. We should point out that 
for the processes with a much more complicated color 
structures, due to the cancellation and enhancement 
of different color structures of the heavy-quarkonium, 
those color states other than the color-singlet state may 
also give sizable contributions. Two such examples for 
the direct hadronic pro duction of Be and Sec can be 
found in Refs. [Tsl. O l36l. [37l| . where the color-octet (cb)- 
quarkonium and the color-sextuplet (cc)-diquark can pro- 
vide sizable contributions up to ^ 10% — 20% to the final 
meson/baryon production cross section. 

For the charmonium production channel — >■ 
(cc)[n] + cs, its total decay width for all the P-wave 
states is 165 KeV, which is comparable to that of t]c or 
J/ip, le. it is about 95% (92%) of that of rj^ (J/ip)- For 
the (c6)-quarkonium production, the total decay width 
for all the P-wave states is about 28% (33%) of that 
of Be {B*) for W+ icb)[n] + bs; and is about 69% 
(47%) of that of Be (B*) for W+ {cb)[n] + cc. Note 
even though — ^ (c&)[«] + cc is CKM suppressed to 
W+ {cb)[n] + bs by iKblVlK^P - 0.2%, it is en- 
hanced by the phase space, since it is easier to generate 
a (cc)-pair than a (66)-pair. So as a combined result, 
the decay width of — )• (c6)[n] + cc is smaller than 
that of W+ -)■ {cb)[n] + bs by only 1/10. For the bot- 
tomonium production, the total decay width for all the 
P-wave states is about 13% of that of ?7b or T. Sizable de- 
cay width for the P-wave quarkonium states shows that 
one needs to take the P-wave states into consideration 
for a sound estimation, especially for the channels of the 
charmonium and the (c6)-quarkonium. 

To show the relative importance among different Fock 
states more clearly, we present the differential distribu- 
tions dT/dsi and dT/ds2 in Figs. (|3l4l5l6p . and the dif- 
ferential distributions dT /dcosOi^ and dT /dcos623 in 
Figs. (|7l8l9llOI) . Here two invariant variables, si = {qi + 
q^)^ and S2 = {qi +92)^, ^13 stands for the angle between 
qi and q^, and 623 stands for the angle between q2 and 
(fa in the iy+-rest frame. The curves for \{QQ')i[^Sq]), 
|(QQ')iP'5i]), \iQQ')i['Pi]) and liQQ'WPj]) are pre- 
sented. Since the difference between the color-singlet 5- 
wave states and the color-octet S'-wave states is an overall 
color factor, the shapes of their curves are the same, so 
we do not present the curves for the color-octet ones in 
these figures. 
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FIG. 3. Differential decay widths dV/dsi (Up) and dT/ds2 
(Down) for — (cc)[n]+cs, where the dashed line, the solid 
line, the diamond line, the crossed line, the dash-dot line and 
the dotted line are for \{cc)i[^Sq]), | (cc)i pSi]) , \{cc)i[^ Pi]) , 
|(cc)ipPo]), |(cc)ipPi]) and |(cc)ipP2l) respectively. 



10"^ 




10"'" I ' ' — ' ' — ' I 

10 10 10 

s,(GeV-=) 



10'^ 




lo'""! ' ' — ' ' — I 

10 10 10 

SjCGeV-^) 

FIG. 4. Differential decay widths dF/dsi (Up) and dr/ds2 
(Down) for — >■ (c&)[n]-|-6s, where the dashed line, the solid 
line, the diamond line, the crossed line, the dash-dot line and 
the dotted line are for | (cfe)i [^So]), | (cfe)i pSi]) , l(c6)i[iPi]), 
|(c5)ipPo]), |(c6)ipPi]) and |(c6)ipP2]) respectively. 



As for the decay channels {cc)[n] + cs and 

— >■ (c6)[n] -I- bs, because s is much lighter than the 
heavy quarks, the largest dT /dcosOi^ is achieved when 
the quarkonium and s move back to back (^13 = 180°) 
or the quarkonium and the outgoing heavy quark move 
in the same direction {623 = 0°), which is shown in 
Figs. ((7I5|) . While for the_decay channels W+ -> {cb)[n] + 
cc and W+ -> (bb) + cb, as shown in Figs. (|9ll0p . the 
largest dF/dcos^is is achieved when ^13 = 0° and the 
largest dT/d 003623 is achieved when 623 = 180°. This 
shows that the maximum differential decay width is ob- 
tained when the quarkonium and the outgoing antiquark 
move in the same direction or the quarkonium and the 
outgoing quark move back to back in the rest frame. 

Considering that the LHC runs at the centcr-of-mass 
energy \/S = 14 TeV with luminosity lO'^^cm^^s^^, 
one expects that about 3.07 x 10^° M^+-events per year 
can be generated Then we can estimate the heavy- 
quarkonium events generated through decays, i.e. 
2.57 X 106 r]c, 2.65 x 10*^ J/^ and 2.40 x 10^ P-wave 
charmonium events per year can be generated; 1.01 x 10^ 
Be, 9.11 x 10'' B* and 3.16 x 10'' P-wave (c6)-quarkonium 
events per year can be generated; 2.74 x 10^ 77b, 2.59 x 10^ 
T and 35 P-wave bottomonium events per year can be 
generated. A small event numbers for the bottomonium 



shows that it is hard to find the bottomonium through 
decays even at LHC. However, it might be pos- 
sible to measure J/tp and Be events through de- 
cays, since W'^ is charged particle and one may identify 
these particles through their cascade decay channels as 
J/t/j M^M^ and Be J/^ + n or Be ^ J/t/j + eVe 
with clear signal. Bearing the situation pointed out here 
and the possible upgrade for the LHC (SLHC, DLHC, 
etc. dll) in mind, the possibility to study the charmo- 
nium and the (c5)-quarkonium via the W boson decays 
is worth thinking seriously about. 



B. Uncertainty analysis 

In this subsection, we discuss the uncertainties for 
the charmonium and the (c6)-quarkonium production 
through decays. 

For the present leading-order calculation, their main 
uncertainty sources include tlie non-perturbativc bound- 
state matrix elements, CKM matrix elements, the renor- 
malization scale and the quark masses m?,, m.e and 
iris. In the literature, the Brodsky-Lepage-Mackenzie 
(BLM) method 1391 or the principle of maximum con- 
formality (PMC) [40] provides a feasible way to derive a 
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FIG. 5. Differential decay widths dV/dsi (Up) and dT/ds2 
(Down) for — {cb)[n]+cc, where the dashed line, the solid 
line, the diamond line, the crossed line, the dash-dot line and 
the dotted line are for \{cb)i[^ So]), | (c6)i pSi]) , \{cb)i[^ Pi]) , 
|(c5)ipPo]), |(c5)ipPi]) and |(c6)ipP2]) respectively. 



precise QCD predictions. The main idea of BLM/PMC 
is to sum all the non-conformal {Pi} terms in the per- 
turbative expansion into the running coupling, and then 
the remaining terms are identical to that of a conformal 
theory and are renormalization-scheme independent. At 
the present, the bound state and CKM matrix elements 
and as emerge as overall factors and their uncertainties 
can be conveniently discussed when we know their val- 
ues well, so we shall not discuss their uncertainties in 
the present paper. In the following, we shall concentrate 
our attention on the uncertainties caused by m^, rric and 
nis, whose values are taken as mi, = 4.90 ± 0.40 GeV, 
rric = 1.35 ± 0.25 GeV and = 0.105 ± 0.025 GeV. 
And for clarity, when discussing the uncertainty caused 
by one parameter, the other parameters are fixed to be 
their center values. 

Typical uncertainties for rric and rrif, are presented in 
TABs. (jVIVI|) . where — )• (cc)[n] stands for the process 
— >■ (cc) [n] + cs, W'^ (c6) [n] stands for the pro- 
cesses — >■ (c6) [n] + bs and W'^ (c6) [n] + cc respec- 
tively. TABs. ljVIVip show that sizable uncertainties can 
be found for varying nib and rric. The decay width will 
decrease with the increment of mf, and rric-, and such ten- 
dency slow down with a heavier quark mass. Taking the 
process W^+(fc) {cc)[n]{q3) + c{q2)s{qi) as an explicit 
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FIG. 6. Differential decay widths dF/dsi (Up) and dT/ds2 
(Down) for — >■ {bb)\n\+cb, where the dashed line, the solid 
line, the diamond line, the crossed line, the dash-dot line and 
the dotted line are for |(fe&)ip5o]), i(6fe)ipSi]), l(b6)i[^Pi]), 
|(65)ipPo]), |(66)ipPi]) and \{bb)^[' P2]) respectively. 



mc(GeV) 


1.10 


1.35 


1.60 


r|(c.),iis„i)(KeV) 


326.0 


174.8 


104.0 


r|{cc)ii3Sii>(KeV) 


336.6 


180.6 


107.5 


r|(cc)i[P-iiiai;e]) (KcV) 


246.9 


165.4 


54.26 


r|(.5),(is„i>(KeV) 


6.33 


6.32 


6.30 


r|(c5)ii3Sii>(KeV) 


5.30 


5.38 


5.46 


^\{cb)-L\P-wave]) (KcV) 


2.48 


1.77 


1.35 



TABLE V. Uncertainties for the decay width of the processes 
— >■ (cc)[n] and — >■ (c6)[n], where \ {cc)i[P —wave]) and 
|(c6)i[P — wave]) stands for the sum of the four color-singlet 
P-wave states for the (cc)- and (c6)- quarkonium accordingly. 



example. One may observe that even though its phase- 
space is slightly affected by rric, i-c. the maximum value 
of si = [qi + q'iY (s7°^ = m^(l - rric/mwY ~ m^) 
remains almost unchanged, the total decay width shall 
be decreased by about 3 — 5 times for various (cc)- 
quarkonium states by varying rric from 1.10 GeV to 1.60 
GeV. Such a big uncertainty is mainly caused by the 
fact that it is harder for an intermediate hard gluon to 
generate a heavier (cc)-pair. More explicitly, at the spe- 
cific momentum region with q^ ~ 4m^ (where qg stands 
for the intermediate gluon momentum), which gives the 
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FIG. 7. Differential decay widths dV/d cos 613 (Up) and 
dr/dcos^23 (Down) for — >■ (cc)[n] + cs, where 
the dashed hne, the sohd hne, the diamond hne, the 
crossed hne, the dash-dot hne and the dotted hne are for 
\{cch['So]), |(cc)ip5i]), |(cc)i[iPi]), l(cc)ipPo]), i(cc)ipA]) 
and |(cc)ipP2]) respectively. 



nib (GeV) 


4.50 


4.90 


5.30 


r|(c5),[iSnl>(KeV) 


8.31 


6.32 


4.89 


r|(cb)i[3Sji)(KeV) 


7.14 


5.38 


4.14 


r|(cb)i[P-ti.ai.el>(KeV) 


2.41 


1.77 


1.33 



TABLE VI. Uncertainties for the decay width of the pro- 
cess — (c6)[n] with varying m^, where |(c&)i[P — wave]) 
stands for the sum of the four color-singiet P-wave states. 



main contribution to the decay width, there is a strong 
suppression factor of (1.60^/1.10''') ~ 4.5 for the S'-wave 
production by varying the c-quark masses from 1.10 to 
1.60 GeV 0. Moreover, one may observe that the de- 
cay widths for P-wave states are more sensitive to the 
quark masses than the case of S'-wave states. Vary- 
ing rus G [0.080,0.130], one may observe that the de- 
cay width of — (cc)[n] -I- cs is almost unchanged 
for the S'-wave states, but shall cause sizable changes 



^ As for the P-wave cases, because of the derivation of the ampli- 
tude over the bound-state relative momentum, this suppression 
factor shall become even bigger. 




FIG. 8. Differential decay widths dT /dcosOia (Up) and 
dF/d cos 6*23 (Down) for — > (c6)[n] -I- bs, where 
the dashed line, the solid line, the diamond line, the 
crossed line, the dash-dot line and the dotted line are for 
m,[\So]), |(c6)ipSi]), |(c6)i[iPi]), l(c5)ipPo]), |(c6)ipPi]) 
and |(c6)i[^P2]) respectively. 

for the P-wave states, i.e. Tii^cc)i[p-wave])\m,=o.o80GcV = 
139.2 KcV, Tncc)i[p-wave]) |m,=o.i05GcV = 165.4 KcV and 

^\{cc)i[P~wave])\ms=0.130GcV = 185.3 KcV. 

Adding all the uncertainties caused by the constituent 
quark masses in quadrature, for W'^ — > (cc)[7i] + cs, we 
obtain 

r|(cc-),[iSo]> = 174.81^^^8^ KeV, 
r|(cc-).[3s,]> = ISO.Gt'rT,' KeV, 

^\(cc)i[P-wave]) = 165.4^112^9 KcV, 
^\{cc)s[S-wave]g) = 44.5l?|;5 >< KcV. 

For + {cb) [n] + bs, we obtain 

r|(.6),[i5ol)=6.32ll:f3KeV, 
r|(c6)ipsi]) = 5.38li;24 KeV, 

^\{cb)i[P-wave]) ~ ^■'^'^^^Tl KcV, 
r|(e6)s[S-^aH9> = ^'^^-Tl >< KcV. 

And for — > (c6) [n] -f- cc, we obtain 

r|(c6),rsol> = 0.546t°;^^«KeV, 
r|(c6)ipsi]> = 0-810to;354 KeV, 
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FIG. 9. Differential decay widths dV/d cos 613 (Up) and 
dV/d cos 023 (Down) for — > (cS)N + cc, where 
the dashed hne, the sohd line, the diamond line, the 
crossed line, the dash-dot line and the dotted line are for 
\{cbU\So]), I(c6)ip5i]), |(cfo)i[iPi]), i(c5)ipPo]), |(c5)ipA]) 
and |(c6)i[^P2]) respectively. 

^\{cb)i[P-wave]) 0-379^o;i56 KcV, 

^\ichMs-u.aue]a) = 0.170t°;J«7° X KeV. 

If assuming the higher excited heavy-quarkonium 
states decay to the ground color-singlet and spin- 
singlet state with 100% efficiency via electromagnetic 
or hadronic interactions, then we obtain the total decay 
width of decay channels, 







iSo]+cS) 


= 524.8l^^^ j KeV 


(33) 


^{W+- 


Hcb)^ 


iSo]+cg) 


= r3.5ttll KeV 


(34) 


^{W+- 


■>(c6)i 


^So]+cc} 


= 1.741 J;^^ KeV 


(35) 






^So]+cb) 




(36) 



where G [0.10,0.30] is adopted. 

Finally, to show how the decay width depends on the 
corresponding quark mass more clearly, we take the chan- 
nel — ?> (cc) [n] + cs as an explicit example, whose dif- 
ferential decay widths with several typical nic are drawn 
in Figs. (|llll2p . The contributions from the color-singlet 
S-wave and P- wave states have been summed up for 
convenience. And in Figs. (|llll2p . the center solid lines 
are for rric — 1.35GeV, the dashed lines are for nic — 
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FIG. 10. Differential decay widths dF/d cos 613 (Up) 
and dT/d cos 623 (Down) for — s> {bb)[n] + cb, where 
the dashed line, the solid line, the diamond line, the 
crossed line, the dash-dot line and the dotted line are for 
\{bbU\So]), |(b6)ipSi]), \ibb)^['P^]), m,['Po]), !(fe6)ipPi]) 
and |(&6)i[^P2]) respectively. 



1.25GeV and the dotted lines are for rric = 1.45GeV re- 
spectively. Similar to the observation from the TabfVl 
large uncertainties are caused for rric varying within the 
region of [1.25, 1.45] GeV. At the specific momentum re- 
gions with Qg ~ 4to^, we have si — m^, which explains 
why there is a peak near si ^ for the differential 
decay width dT/dsi as shown by Fig.pTj). 



IV. CONCLUSIONS 

In the present paper, we have made a detailed study 
on the heavy-quarkonium production through W bo- 
son semi-inclusive decays, — > (cc)[n] + cs, — !> 
(c6)[n] -I- cs, {d>)[n] + cc and W+ -> {bb)[n] + cb 

within the NRQCD framework. Results for eight quarko- 
nium Fock states, i.e. |(Qg')i.8['^o]), \iQQ')i,s?Si]), 
|(QQ')iP-Pi]> and \iQQ')i[^Pj]) have been presented. 
And to provide the analytical expressions as simply and 
compactly as possible, we have adopted the improved 
trace technology to derive the Lorentz invariant expres- 
sions for the W boson decay processes at the amplitude 
level. Such a calculation technology shall be very helpful 
to derive simpler results for the processes with massive 
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FIG. 11. Uncertainties of dV/dsi (Up) and dT'/ds2 (Down) 
for — >■ (cc) + cs, where contributions from the color- 
singlet S-wave and P- wave states have been summed up. 
The dashed line, the solid line, the dotted line are for rUc = 
1.25GeV, 1.35GeV and 1.45GeV respectively. 



spinors. 

Numerical results show that the P-wave states in ad- 
dition to the S-wave states can also provide sizable con- 
tributions to the heavy-quarkonium production, so one 
needs to take the P-wave states into consideration for 
a sound estimation. More explicitly, for the charmo- 
nium production channel — > (cc)[n] -f cs, the to- 
tal decay width for all the P-wave states is 165 KeV, 
which is comparable to that of 77c or J/il). For the (c6)- 
quarkonium production, the total decay width for all the 
P-wave states is about 28% (33%) of that of (B*) for 
W+ {cb)[n] + bs; and is about 69% (47%) of that of 
Be (B*) for W+ {cb)[n] + cc. If all the low-lying ex- 
cited states decay to the ground state \{QQ')i[^ So]) with 
100% efficiency, we can obtain the total decay width for 
the {QQ') production through W'^ decays as shown by 
Eas. (|33l34l35l36p . At the LHC, due to its high colhsion 
energy and high luminosity, sizable heavy-quarkonium 
events can be produced through W'^ boson decays, i.e. 
2.57 X 10^ 77c, 2.65 x 10^ J/'f and 2.40 x 10^ P-wave char- 
monium events per year can be obtained; and 1.01 x 10^ 
Be, 9.11 X 10'' Bl and 3.16 x 10'' P-wave (c6)-quarkonium 
events per year can be obtained. So these channels may 
be an important supplement for other measurements at 
the LHC. And large event numbers for the higher excited 
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FIG. 12. Uncertainties of dr/dcos6'i3 (Up) and dV / dcos92i 
(Down) for — >■ (cc) [n] + cs, where contributions from the 
color-singlet S-wave and P- wave states have been summed 
up. The dashed line, the solid line, the dotted line are for 
771c = 1.25GeV, 1.35GeV and 1.45GeV respectively. 



states may cause themselves to be deserving of careful 
study. 
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APPENDIX: AMPLITUDES FOR 

W+{k) (cQ')N(g3) + Q'(g2)%i) 

We shall only list the results for the first type of semi- 
inclusive process W'^{k) — >• (cQONl'Za) + Q'{Q2)s{qi)- 
The results for the second type of process W'^{k) ^• 
{Qb)[n]{q3) + c{q2)Q{qi) is similar, 
o Short notations. 

To shorten the notation, we first define some dinicn- 
sionless parameters 

777.Q' 771, M 

n = , 7-2 = — -, ra = , 

mw mw mw 
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X = z ^ ,y = 



For \{cQ')i[ Sq]) , there are three Lorentz structures, 



2mw' 



mw +rns -S3 S3 ~ mni - 

z = z ,1* — 



2mw^ 
si - - 
2mw^ 



2mw^ 



S2 - mq,^ - nis'^ 



2mw^ 



, (37) 



93 ■ £(fc) „ g2 ■ £(fc) 

-t>l — , iJ2 — : 

1 

B3 = T£(fc,'73,'72,e(fc)), 



(40) 



where si = (gi + ga)^, ^2 = (gi+92)^ and sg = (92+q3)^ where e{k,q3,q2,e{k)) = e{a, /S, p,a)kaq3i3q2pea{k). 
which satisfy the relation si + S2 + S3 = jtivj/^ + 



mQi^ + Af ^. As for the normahzation constants, we have Lorentz structures 



For KcgOipS*!]) and |(cQ')i[^A]), there are eleven 



Li,2 = l/{2'mw\/w =F rir2) and 



™vF'^'\/?'i^''2^ — — r2^y'^ ~ ri^z^ + 2wyz 

The short notations for the denominators are 
1 1 



dii 
d2 
d2i 

d22 



mq''^ (92 + 532) 



,0 \2 ' 



niw 



sl-mQ.^{q2 + ql2Y' 



1 



1 



('Zl+g2+93"2)'-m.2 (g2+g3°2)'' 

■2 1 



((gi + g2 + 9!^2)'-™.')' fe + 93°2)'' 



((gi + g2 + g3°2)'-™.')' (92 + Q3°2)^' 



(38) 



where q^^ = ff^s and = ^93- Finally, the follow- 



ing relations are useful for further simplification, 

X + y + z = 1, u. + i) + w~l, 

u + w + rg = .X, w + u + r\ = w + v + = 

o Basic Lorentz-invariant structures. 



.(39) 



J 



Bi = £(93) • e(fc), B2 = ^-£(fc, 93, e(fc), £(93)), 



B3 = 



Bi = 



-e(fc,(72,e(fc),e(q'3)), 
-e(93,'72,e(fc),e(g3)), 



_ k- e{q3)q3 ■ e(fc) fc • £(93)92 • e(fc) 

^5 — 9 , — 9 • 



o 92 ■ £(93)93 ■ e(fc) o _ 92 • e(fc)92 • £(93) 

ii? — n , J^R — n 



m 



W 



m 



w 



B. = 



B 



10 



-e(fc,93,92,e(fc))fc • £(93), 
-£(fc,93,92,£(93))93 ■ £(fc), 



•■w 



Bii = -^e{k, 93, 92, £(fc))92 • £(93) 



(41) 



where 2(93) stands for the polarization vector relating to 
the spin or the orbit angular momentum of |(c(3')ipS'i]) 
and |(c(5')i[-^Pi]) respectively. 

For \{cQ')i[^Pj]) with J = 0, 1, 2, there are thirty-four 
Lorentz structures. 



Bi 
B5 
Bs 
Bii 

BlA 

Bn 
B20 
B23 



^—92 • e{k)ei^, B2 = ^—93 • £(fc)£Qa7 B3 = ^— 92a£/3(fc)ea/3: B4 = ^— fcQ£/3(fc)e^ 
mw niw mw ^ mw 

J 

e(fc,93,£(fc),a)92/3 



/J' 



^^£(fc,92,93,£(fc)), Be = —^e(fc, 92,93, a)£/3(fc), S7 = — ^ 

m'^ 



"'^ e{k,q2,a,l3)q2 ■ e{k), Bg = — |^e(/c, 93, £(fc), a)fc^, Bio = — f^£(fc, 92, £(fc), a) 



"'^£(fc,93,a, /3)93 • £(fc), B12 = — |^£(fc,92,a, /3)93 • £(fc), B13 = — 



e{k,q3,a,l3)q2 ■ e{k), 



-£(fc,92,£(fc),a)92/3, -B15 



£(92,93,£(fc),a)92/J, B16 



£(92,93,a,/3)92 • £(fc), 



«£ 



Q/3 



e(92,93,a,/3)93 ' £(fc), -Big 



l£ 



(92,93,£(fc),a)fc/J, -B19 



i£ ' 



a/3 



£(fc, £(fc), a, (i), 



le 



a/3 



£(92,£(fc),a, /3), B2 



£(93, £(fc), a, /3), B22 = -^e{k, 92, 93, a)fc/393 • £(fc), 

mw m 



w 



—^£(fc, 92,93, Q;)92^92 ■£(fc), B24 = — ?^£(fc, 92, 93, a)fc/392 ' £(fc), -B25 = 
TO^ TO^ 



92,93, a)92/393 • £(fc), 



13 

XS'^ IS"^ %S'^ 

B26 = —T^£(fc, 92,93, e(fc))'?2a<?2/3, S27 = — ^e(fc, (J2 , 93 , e(^) fc/3 , ^28 = — ^^(fc, 92 , 93 , e(^))^a92/3 , 
771^ TO^ TO^ 

i?29 = -^kakpqs ■ e(fc), = -^^092/393 • e(fc), -831 = -^kakisq2 ■ e(fc), 

B32 = -^fca92/392 • e(fc), i?33 = -^92a92/393 ' S34 = -^92q92^92 ' e(fc)- (42) 



o Non-zero coefficients for \{cQ')i[^So]) . Non-zero coefficients and are 



or 7/2 

A} = ' ^ (d2((2r3 ~ 4r2)ri2 + (r32 - 4r22 + 2?/ + 2z)ri - 2(r2 - r3)(M - y)) 



+2di(~2ri^ - 2(r2 - ra)^' + {r^' + 2r2r3 - 2u + 2y)ri + Ar^u - 2r^y)), (43) 
A\ = ^^'"^^^' (rf2(ri + r2)(r32 + 2r2r3 - 2x) - 2di(-2r33 + r2r32 + 2a;r3 + 2r2U + r^{r^'' - 2v))), (44) 



^1 ^ + _ 2di(ri - rg)), (45) 

\A3 



= - ° (-4d2?/^ + (2di(-2ri2 + 2(r2 + r3)ri + r32 - 2r2r3 + 2u) + rf2((2ri - 2r2 + r3)(2r2 + r3) 

+4ii + 2))j/ + 2(2diri(ri - r3)(l - 2x) - d2{u + ri(ri - r2 - r3)(x - 1)))), (46) 

= (d2((x + 2y - 3)r32 - 2r2(x- + 2y - l)r3 + 2x - Axy) 

-2di((a; + 2y- \)rz^ + 2(ri + r2)a-r3 - 2u - 2ri(ri + r2)x + 2u{x + 2y))), (47) 
A% = !^Jl!^!^(4di(ri + r2)(ri - r^) + 2^1(73^ + 2u) - d2{r^^ + 2r2r3 + 4y - 2)). (48) 



Non-zero coefficients A| and A| are 



A{ = - ^^'"^'^' (rf2(2(2r2 + r3)ri2 + (73^ - 4r22 + 2u + 2z)n + 2{r2 + r3)(u - y)) 



+2di{-2ri' + 2(r2 + rj,)ri^ + (r3^ - 2r2r3 - 2u + 2y)ri + 4r3U - 2r3y)), (49) 
Al = ^^'"!!!^^^' (2di(2r33 r2r32 - 2a:;r3 + 2r2W - r^{r:i'' - 2v)) - rf2(ri - r2)(-r32 + 2^273 -|- 2x)), (50) 

\A3 



fr^ 



Al^ ' ^ (d2(r2-ri) + 2di(ri-r3)), (51) 



Ai = !^^!^!^(-4diri(ri - r3)(2a; - 1) - 2d^{-r^ - 2(ri + r2)r3 + 2ri(ri + r2) - 2u)y 

+d2(-4y2 + (^3^ + 2rir3 - 4r2(ri + ra) + Au)y + 2y~2u- 2ri(ri + r2 - r3)(x - 1))), (52) 

Ai = - '^^''^^\ d2{{x + 2y- iW + 2r2(y - z)r3 + 2a; - 4xy) 

-2di((y - z)r32 + 2(ri - r2)xr3 + 2(ri(r2 - ri)x + u{y - z)))), (53) 

4' = "'^'^'i!° (d2(r3^ - 2r2r3 + 4j/ - 2) - 2di(2ri2 - 2(r2 + r^^^ + r^"" + 2r2r^ + 2u)). (54) 
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o Non-zero coefficients for \{cQ')i['^Si]) . Non-zero coefficients Aj and Aj are 



Al = / (ti^(2(r-32 _ 2r2r3 - 2x)ri^ + {ra^ + (u - 4r2^ -x + 2z)r3 - 2r2(M + x))ri - 2r2^ 



+2(r3^ - 2x)iu - y) + r2r3(2y - 3u)) - 2d^r^{r2U ~ r^v)), (55) 
Al = ^^'"!^'^' (d2(ri(4ri + 2r2 + r3) + 8(u - y)) - 2d,r^r^), (56) 



= + _ 2r,)r, + 2x), (57) 



2L m 

= ^—^ (2di(ri + r2)r3 + d2((ri + r2)(2r2 + rs) - 4x + 4)), (58) 

\/r3 



2L m 

= (d2(ri(4ri + 2r2 + r3) + 4(u - y)) ~ 2dirir3), (59) 



Al = (^^^ ^1 = 4L,d2mw'^'V^{ri + r^), (60) 



Al=A\, = -'-h!^^, (62) 



Af = - " ^ (-2dir3(M + ri(ri + r2)x) + 2dir3(r32 + 2u)y + d2((r-3(x - 2) + 2r2.T)ri2 + (2a;r22 + r3(-3a; 

-4y + 2)r2 + 2(r32 - 2x)z)ri + 4r2V3y + 2r2(it - 2{u + x)y) - r3(2yu + u+ {r:^'^ + 2)y - 4y(x + y)))), (63) 
= - ^°'"^!^' (2dir3(ri(ri + r2) + 2if) + rf2((r3 - 2r2W + ('^2(^3 - 2r2) - 4x + 4)ri - 2(2r2 - r^){u ~ 2y)p,4) 

At = ^°"'^!^' (2dir3^ + d2(r33 - 2r2r32 + (4y - 2)r3 + 4(ri + r2)x-)), (65) 

= - ^"'''^!^' (-2rfir3 + 4di(x + y)r3 + d2(4ri + (2r2 - r3)(l + 2z))), (66) 



/y ™ 9/2 

^i' = ° ((2dir3(ri(ri + r2) + 2w)) + (d2((r3 - 2r2)rf^ + (r2(r3 - 2r2) + 4z)ri + 2r3W + 4r2y - 4r3y))|p7) 



^3' ^ ^!2!^!^^^(2dir32 + d2(r32 + 2r2r3 + 4y - 2)), (68) 
= ^°"'^!^' (2dir3(2x + 4y - 1) + d2(r2(2 - Ay) - r3(l + 2z))), (69) 



, = 2^^W^ ^3' ^ „ 2d2W/- ^3; ^ (^^\ ^ _2d2iW^_ 

^2 ^2^^ \dir3J L2y^ 



r 



Non-zero coefficients Aj and are 



9 m ^/ 

^2 ^ ^ ^ {2dir3{r2U + nv) + d2{2{r3^ + 2r2r3 - 2x)ri'^ + {r^i^ + {u ~ 4r2^ ~x + 2z)r3 



-|-2r2(M -t- x))ri - u((r2 - 2r3)(2r2 + r^i) + Ax) - 2(r3(r2 + r^i) - 2x)y)), (71) 



2 2L2TO^v^^^ 



= ^ (d2(r-i(4ri - 2r2 + r3) + 8(u - y)) - 2dirir3), (72) 
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2 Ad2L2mw'^^'^ , 



Ai= VJ'^ (rir3-(r2 + 2r3)r3 + 2x), (73) 



2 2L2Tnw'^^'^ , 



Ai = (2di(r2 - ri)r3 + d2((ri - r2)(2r2 - rg) + 4(x - 1))), (74) 



2 2L2"lH'^/^, 



= ^ (d2(r-i(4ri - 2r2 + r^) + A{u - y)) - 2dinr3), (75) 



= (^||) - 4d2i2mH.^/2(ri - r2)V^, (76) 
^2 ^ ^^™'^' (2di(3ri - r2)r3 - ^2(^1 - r2)(2r2 - r3)), (77) 



Af = (2dir3(yr3' + ri(r2 - ri)x + u(2y - 1)) + d2((r3(a; - 2) - 2r2a:)ri2 + (2a;r2' + r3(3x + 4y - 2)r2 

+2(r32 - 2a;)z)ri - 2r2U + 4r2V3y + 4r2(u. + x)y - r3{2yu + u + {rs^ + 2)y - Ay{x + y)))), (79) 
= ^^p^{2dir3{n^ - + 2u) + d2{{2r2 + ny^"" - (2r22 + r3r2 + Ax - 4)ri + 2(2r2 + r^){u - 2y))), (80) 



Ai = - ^"'"^!^' (2rfir3^ + d2(r3' + 2r2r32 + {4y - 2)r^ + 4(ri - T2)x)), (81) 
Ai = ^^p^(-2d,r3 + 4di(x + y)r3 + rf2(4ri - (2r2 + r3)(l + 2^))), (82) 
Ai = ° (2dir3(ri^ - rzri + 2u) + d2((2r2 + r3)n^ - (r2(2r2 + r^) - Az)r^ + 2r3U - 4(r2 + r^)y)), (83) 



Ai = - ^°"'!;f' ^ (2dir3^ + d2(r3^ - 2r2r3 + 4y - 2)), (84) 
Ai - -^^^^^(2rfir3(2a; + 4?; - 1) - rf2(r3(l + 2z) - r2{4y - 2))), (85) 
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Ai = _Ml^W^V?(^ _ ^4' ^ 2d2A^om^V^ _ r d_2r2\ ^ ^2_d^r2N^^rn^ ^ 



o Non-zero coefficients for \{cQ')i[^ Pi]) . Non-zero coefRcients and A^ are 



2L m 

= \ T/2 (c?2r3'(-2r3ri'^ + (2a; - 2r2r3)ri2 + (r2(2a; - v) + 2r3(y - u))ri + r2^u + 2ux - 2xy) 

-2d2iri2((ri2 - r2^ + l)r3 - 2rix)(2r3ri2 + (2r2r3 - v)ri + r2W + 2r3U - 2r3y)), (87) 
A\ = -'^h^^^{2d2i{{r^ - r2^ + l)r3 - 2rixW + d2r^{2r^ + r2ri + 4u - Ay)), (88) 

.1 4d2Limvy5/2y^ 

^3 = 5 — (?-3 -a:, (89 

A\ = '^^'Y'X''\ d2r^{r2{ri + r2) - 2x + 2) - 2d2iri2(ri + r2){{r^ - r2^ + l)r3 - 2rix)), (90) 
^5 = -^^^?^(2d2i(3r3ri3 + 2(r3^ + 2^. - x)^^ + (-3r3r22 - 2r32r2 + 4z.r2 

^1 '3 

+f3 + 4r3U - 6r3?;)ri + 2r2rz{y ~ 2u))rx + d2rz^(r2rx ~ 2w)), (91) 



■;= (ri + r2)((ri + r2)r3 - a;), (92) 

JlUmw^l ^ r2)((ri2 - r^^ + l)r3 - 2ria-)ri2 + 4d22(2r3ri3 - (rg' + 2vW 

+(-2r3r2^ - rg^ra + 2ur2 + rg + 2r3U - 3r3y)ri + r2r3(y - 2M))ri2 + r3(4(iii(-2ri3 - 2(r2 - r-i)r^^ 

+(r32 + 2r2r3 - 2ii + 2y)r^ + 4r3W - 2r3y)ri2 + ^2^2(^1 + r2)r3)), (93) 

' ^ -id22iri+r2)iiri+r2)r3~x)-dui-2r-/ +r2r3^ + 2xr3 + 2r2U + n{r3^ -2v))), (94) 



/f3 

^hlp^{d2r3-2d2in'in+r2)), (95) 
M2Limw^^^y^ 

2 — ' (96) 



(d22(r-i+r2)-2dii(ri-r3)), (97) 



/r3 

2^2^1^33/2" 

+d2r3^iir2X + 2r3z)ri^ + x{r2^ + 2x - 2)ri + 2(r2r3 + x)yri + r2{u - (rg^ + 2{u + x))y))), (98) 
iVomvy'^/^ , , , , N , o_ 2 2 , IN o \ 2 



(2d2i((ri2 - r2^ + l)r3 - 2rix){2rir3 + u + nin + r2 ~ 2r3)x - (r3(2ri - 2r2 + rs) + 2u)y)r{' 
- 2r3z)ri^ + a;(r2^ + 2x - 2)ri + 2{r2r3 + x)yri + r2 
(2d2i(ri(ri + ra) + 2u){{ri'^ - r2^ + l)r3 - 2ria;)ri^ 



2L2ri2r33/2' 

+d2r3^(r2ri2 + (rs^ + 2x - 2)ri + 2r2(M - 2y))), (99) 

^°'rr'^'2^ (2^2i((ri" - 1-2' + l)r3 - 2r^xW + d2r2r3' - 2^2(^1 + r2)x), (100) 

_ jVomvyV2 2 _ ^^2 ^ ^^^.^ _ 2rix)(l - 2z)ri2 - d2r3'(r2(l + 2z) + 2n)), (101) 

2-/j2f 1 ^3 ' 

^ r(2rf2i'^i'(4r3y' - 2(r3(2ri2 + + r2(r3 - 2r2) + 1) + Ar,u)y 



2i2ri2r33/2" 

+ri(4u + (ri + r2)(r3(ri2 - r2^ - 4x + 3) + 2rix))) - d2r3^{r2ri^ + {r2^ - 2z)ri - 2r2y)), (102) 

jvr ™ 7/2 

° (2d2iri2(ri2r32 - (r3 + r2(r2r3 + 2x - 2))r3 + 2ri(2y - l)r3 + 2x - 4(r2r3 + x)y) - d2r2r33)(103) 



2£2ri2V^ 

_ ^Q-Tfiw (-Sd22xri^ - 4:d22r3ri^ - 8d22r2xri'^ + 8d22r3xri'^ - 4d22?'2?'37'i^ - 8d22uri^ 
2X2^1 r3'''^ 

+8d22r2r3xri^ - 8d22r3y^n^ + 2d2i((ri^ - + l)r3 - 2rix){2y - l)?'!^ + 4diir3{{-2ri^ + 2{r2 

+r3)ri + r32 - 2r2r3 + 2u)j/ - 2ri(ri ~ r3)(2a; - l)W + ^27-2^3' + 2(2d22ri2(r3(2ri2 - 2r22 

+ (ri + r2)r3 + 1) + Arm) - d2r2r3^)y), (104) 

2iVoTOH/ ''/2 , , , o , , , , „ , 

P — ^^(022(^-3 -ri{y- z)r3 + r2{y ~ zjrs ~ x + 2xy) 

+dn{{y - z)r3^ + 2(ri + r2)xr3 - 2m - 2ri(ri + r2)a: + 2u{x + 2y))), (105) 



„ ^(2d2iri ((ri+r2)r3 + 2?;-l)-rf2(ri+r2)r3), (106) 

/2 

(2d2i((ri2 - r22 + l)r3 - 2rix)n^ + d2r2r3^), (107) 



L2ri'^y/r^ 
L2ri^r3^/'^ 

-(rfii(r-3^ - 2(ri + r2)r3 + 2ri(ri + ra) + 2w) - cf22((ri + r2)r3 + 2y - 1)). (108) 



Non-zero coefficients and are 



17 

\ T/2 (2^21 ((ri' " + l)r3 - 2ria;)(2r3ri2 - {v + 2r2r3)n - + 2r3U - 2r3y)n^ 

+d2r3^(2r3ri^ - 2(r2r3 + a;)ri2 + (r2(r3^ + u + x) + 2r3{u - y))ri - u{r2^ + 2x) + 2xy)), (109) 
2L2mw^^'^ ,^ , ,,2 2 , IX o \ 3. , , 2(r, 2 



7'-j^2^^3/2 



(2d2i((ri' ~ + l)r3 - 2rix)ri' + rf2r3'(2ri' - rzri + 4u - 4y)), (110) 



(r3 -a^), (111) 



2L2mw^/^ 



(d2r3'(r2' - rir2 - 2x + 2) - 2d2in''{n - r2){{ri^ - ^2^ + l)r3 - 2rix)), (112) 
(2d2i(3r3ri3 + 2(r32 + 2u - xYi"^ + (-3r3r22 + 2r32r2 - 4Mr2 + 73 + 4r3M 



2i2mH^5/2 



-6r3y)ri + 2r2r3(2u - y))n'' + d2r3'{2n^ - r2T^ +2u^ 2y)), (113) 

8d2iL2mw^^'^ , \t , , \ fTiA\ 

— {ri~r2)(-rir3 + r2r3 + x), (114) 

2 3/2 (2rf2i(?'i - r-2)((ri2 - r22 + l)r3 - 2nx)n^ + 4d22{2r3n^ + {r^^ + 2u - 2xW + (-2r3r2' 
+r3^r2 - 2ur2 + 7-3 + 2r3M - 3r3y)ri + r2r3(2M - y))r^,^ + r3(4dn(-2ri3 + 2(r2 

+r3)ri2 + (r32 - 2r2r3 -2u + 2y)r^ + 4r3U - 2r3y)r{^ + d2r2(r2 - r^^^)), (115) 

^^'"^'^' (^22(^1 - r2)(rir3 - r2r3 - x) + dii(2r33 + r2r32 - 2a;r3 + 2r2U + r^i-iv^^ -2u + 2x))), (116) 
4L2mvi/^/^ 



-(2d2i(r2-ri)ri^ + d2r3), (117) 

4d22.2mvi/5/2^ 

2 ' ^^^^1 

8L2™w^/' (d22(r2 - ri) + 2dn(ri - r^)), (119) 

(2^21 ((ri^ - r22 + l)r3 - 2r^x){-xr^^ + (r2X - 2r3z)r^ + r3(2r2 + r3)y + u(2y - \y)r^^ 



IT3 
iVomv^V2 



2Liri2r33/2 

+rf2r3^((r2X - 2r3z)r{^ - x(r2^ + 2x - 2)ri + 2(r2r3 - x)yri + r2(w - (r^^ + 2(u + a;))y))), (120) 

o'^°TT/2 (2^21(^1' - + 2")((^i' - ''2' + 1)^3 - 2rix)ri2 
2Liri^r3'^' ^ 

+d2r3'(-r2ri2 + {r2^ + 2x - 2)ri - 2r2(w - 2y))), (121) 

2£iri2 
2Liri2r33/2' 

r (2^21 (47-3 j/'^ - 2(r3(2ri2 + r3ri - r2(2r2 + r^) + 1) + 4riu)y + ri(4u + (n 

-r2)(7^3(ri2 - r2^ - 4x + 3) + 2r^x)))r^^ + d2r3'(ri((ri - r2)r2 - 2x + 2) - 2(ri + r2)y)), (124) 

^ (rf2r2r3^ + 2(i2iri2(r3(r3ri2 - 2ri - r3 - r2(r2r3 - 2x + 2)) + 2x + 4((ri + r2)r3 - a;)y)),(125) 
2X17-1^^/7^ 



(d2(r2(r3' - 2x) + 2rix) - 2d2ir-i'((ri^ - T2' + l)r3 - 2rix)), (122) 
(2d2iri2((ri2 - T2^ + l)r3 - 2ria:)(l - 2z) ~ d2r3\2n ~ r2(l + 2z))), (123) 



-8d22a;7'l - 4d227'3''l + 8d227'2a;7-i + 8d227'3a;7'i + 4(i227'27'3^l" 



2Liri2r33/2 

-8d22U7-i3 - 8d227-2r-3a;ri3 - 8d22r-3y^?^i^ + 2d2i((ri2 - + l)r3 - 2rix)(2y - l)ri2 
+4diir3((r3^ + 2(ri + r2)r3 - 2ri(ri + ra) + 2u)y - 2ri(ri - r3)(2x - l))ri^ 

-^27-27-3^ + 2 (2^22 (7-3 (2ri 2 + r3ri - r2(2r2 + rs) + 1) + 4riu)ri2 + d2r2r3^)y), (126) 
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Aj = '^^o'mw — (■^22(r3^ - ri{y ~ z)r3 - r2{y - z)r3 - x + 2xy) + - z)r3^ 



+2(ri - r.2)xr3 + 2(ri(r2 - n)x + u{y - z)))), (127) 



Af - ^°"'f'/! (2d2i(rir3 - rarg + 2y - l)ri2 + ^^(ra - r^rg), (128) 



,4, A^ofTHv^/' 



^10 = - r 2 ,/2 (2rf2iri"((ri" - + l)r3 - 2rix) - ^2*^2^3"), (129) 
Afi = - ^^°"']l^^' (rfii(2ri^ - 2(r2 + r3)n + rg' + 2r2r3 + 2u) + d22(-nr3 + r2r3 - 2y + 1)). (130) 



(132) 



o Non-zero coefficients for \{cQ')i[^Pj]){J = 0,1,2). A}(|(cg')ipPi])) = for i = (1-4), 

It is noted that £°'^ is the symmetric tensor and e^^ j = (1, 2, 5, 26, 27, 29, 31, 33, 34) 

is the anti-symmetric tensor, and the fact that e^^^ = 
^aa = 0. so the terms involving the foUowing coefficients 
have no contributions to the square of the amphtude, and 

practically, we can safely set the coefficients before them ,i , ^.^ ri „ ^ 

to be zero A^HWU^ P,])) ^ for z = (1-4), 

j = (1, 2, 5, 8, 11, 12, 13, 16, 17, 19, 20, 21). 

^;(|(cQ')iP^o])) = for I = (1-4), 
j = (8, 9, 11, 12, 13, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25) 

(131) Non-zero coefficients A] and for \{cQ')i[^ Pj]) are 



or ^ 5/2 

Al = \ ^_ ((ri + r2)(2(i2i(»-3ri^ - 2a;ri - r22r3 + r3)ri'^ + d2{r3^ - 3rir32 + 2ri2r3 - 2rir2r3 - xr3 + 2rix)) 



+di(8r3ri3 + 2(r32 -|- 4r2r3 -f- 2u)r{^ + (-5r33 + 2r2r32 - 2Mr3 + 4a;r3 + 4r2u)ri - r2r3(r32 ^ 2w))), (134) 

A\ = '^^^^^.T {d2iAr3ri^ + (8r2r3 - 6r32 - 4z;)ri3 - (3r33 + 4r2r32 + 2(-2r22 + u - 3x + y + l)r3 - Ar2u)ri^ 

+^•3(^3^ + (-2r2^ - 2w - 2x + 3y + l)r3 - 2r2(w + y))ri + r2r3^y) - ri{2d2iri{{ri'^ - + l)r3 
-2rix)(ri(2(ri + r2) + rg) + 2(m - y)) + dir3{l2r3ri^ + 2(r3' + 2r2r3 + 2u)ri - r3(r3' - 6u + 4y)))), (135) 

1 SLlniiY^^'^ , , , \//2 2n 

A3 = — [diir3{r2U ~ rivj - d2i(ri + r2)((ri - r2 + l)r3 - 2ria;) 

v''3 



-f(i22(2r3ri3 -f- {-r3^ + 4r2r3 + w + x)ri^ - (r33 + r2r32 + {-2r2^ - u - 2x + y + l)r3 

-r2{2u + x))ri + r2^u + r2r3{u - J/) - 2(r3^ - x){u - y))), (136) 
AX = — 3r3ri-^ + -r3^ + 4r2r3 +u- x)ri^ + r3r2-' + -r3^ + 2u + x)r2 



■'r3{r3^ -'U~2x + y))ri -f r2^u + r2r3('u ~ y) - 2(r3^ - x){u - y)), (137) 

(138) 





2LiTOh/^/%j , , 
5-^^(02(^-1 + r2)(r3 - 


2ri) 


- 4dirir3), 


A\- 


^L^mw^l'^ 2 

= (a22(3ri + r2rx - 


f 4u 


- 4y) - diirirg). 




8(i2i£iTOH/^/%„ 2 , 
— (3ri +r2ri-f- 


4u - 


-4y), 




/ ^21 \ _ 8d2liimvK^^^ 


((ri- 


f r2-2r3)r3-|-2x). 


^11 — 


4^2^1771^1/^/^ 
7^3 







(139) 
(140) 
(141) 
(142) 



2 — ('^i + ^2), 



ir-i 



(4d2i((ri^ - r-s^ + l)r3 - 2rix)ri2 + d2{ri + 3r2)r3^), 
idii{ri + r2)r3 + d22{-ri'^ + r^ri + r2{r2 + r^) - 2x + 2)), 



{d2r3{ri{ri + ra) - {n + 2r2)r3) - 2d2iri'^{{ri'^ - + l)r3 - 2rix)), 



-(-ri^ + ran + r2(r2 + rg) ~ 2x + 2), 



8d2lilTOH'^^^ , 2 



/r3 
2LimH'5/2 



(2d2i(rir3 + 2u)((ri2 - + l)r3 - 2nx)ri' 



+?'3^(''i?'3(rf2?'i + dira) + {2diri + d2{ri + 3r2))it)), 

2LiTOh/^/^ , , „ , , 2 2 , i\ o \ 2 , 2 



(2d2i((ri + r2 - 2r3)r3 + 2x){{ri^ - r2' + l)r3 - 2rix)ri' + r3^(d2?-3(?'i(r-i + r2) 
lr-3) + di(ri + r2)(r32 + 2u) + d2{ri + 3r2)a;)), 
(2d2i((ri2 - ^2^ + l)r3 - 2rix)(ri(2(ri + r2) + r3) + 4w - 2y)ri2 



^1 ^ r3'^/^ 

-(ri + 3r2)r3) + di(ri + r2)(r32 + 2u) + d2{ri + 3r2)a;)), 
2£i?nty5/2 2 „2,1^„ ^o^„ , „ ^ , . ^ . o„a. 2 

+r32(diri(r32 + 2u) + d2{{2r2 + r3)ri^ + (2r2^ - r3^ + x - z)rx + r2(4M - y)))) 
^ ^ ^ ^ ^ ^ 16LimH.5/2 

(^22 ^21 ^22 

8d2iiimH/^/% 2 , o ^ 

-(ri + r2r\ - 2w), 

(rf2i(fi + r2){ri - r2 - 7-3) + diirir3 - ^22(^1^ + r2ri - 2w)), 



lr3 



^ = ^^SL,m^'/'V^3in+r2), 
"22 ^"11 

SZ/irfHy (c;^-^(3^^ + r2)r3 + d22(-»'i^ + ^3^1 + r2(r2 + r3))), 



:(2d2i((ri' - r2' + l)r3 - 2rix)(y - z)ri^ + d2(2r3(y - z)ri' 



2L2ri^^' 

+ ((-3a; - 6y + 5)r3'^ + 2r2{y - z)r3 - 2x + 4xy)ri + r3{{x + 2y- 2)r3'^ + x- 2xy)) + di(4r3(3a; 
+4y - 2)ri2 + 2((y - z)r3^ + 2r2xr3 + 2u{y - z))ri - r3(r32 + 2u){y - z))), 

^ ^omw^/^ 2 _ ^^2 ^ _ 2rix)((r3 - 2r2)y - 2nz)n^ + dir3(4r3(2x + 3y - 1)^^ 

+2(r32 - 2r2r3 + 2u)yri ~ r3{r3^ + 2u)y) + d2{4:Xri^ + {Ar2X - 2r3{hx + 2y- 3))ri^ 

+ ((6x + 2y~ 5)r3'^ - 2r2(a- - 4y + 1)7-3 + 4u - 8My)ri2 + r3{-4y'^ + (-4r2^ - 4r3r2 + 3r3'^ + I2u + 2)y 

+r2r3 - 6ii)ri - r32((-27-22 + r32 - 2y + l)y + w(4y - 2)))), 

^ — ° (d2i((?'i^ - 7'2^ + IjT-s - 2rix){y - z) + diir3{-yr3^ + u + ri{ri + r2)x - 2uy) 
L2y/ri 

+d22{xri' + 7'3(y + 3z)ri^ + {-xr2^ + r3(3y - z)r2 + u - 2{x - l)x - 2{u + x)?; 
-r32(y + 2z))ri - 2r22r3?; + r3(r32 + 2u + 2^; - l)y + r2(-7;7-32 -u + 2{u + x)y))), 

2d2imw''' Nq 4 \ 3 / 2 n \\ 2 / 3 2 

■p= {-r3ri + (3.T - 7'27'3)ri +(2r2X + r3(r2 - u - x + 2z))ri + {r3r2 ~ xr2 

^2^A3 

+r3{x + 4y - 2)r2 + u + 2(u - + 2x - 2{u + x)?/ - 7-32(2; + 2z))ri + r22r3(u - 2y) 
+r2{-yr3^ + 2{u + x)y) + r3{{r3^ + l)y - 2{x + y)y + u{2y - 1))), 



+d2{2r3n^ - (Sra^ + 2r2r3 + 4y - 2)ri + rgCrg^ + 2y - 1))), 

2d2imiy^/^iVo 2 2 , IN o N 
((r-i - r2 + Ijra - 2rix) 



L2y/ri 

-(rfiir3(ri(ri + 7-2) + 2u) - d22{ri^ + (2r2 - r3)ri' 



L2y/rE 

+ {r2{r2 - 7-3) + 2{u + x - l))ri + 2{r2U - r^u - 2r2V + r^y))), 

2d2imw'/^No_^^^3 + (2r2 - r^y^ + ^2' - rsr2 + 2iu + x - + 2(r2{u - 2y) + r,{y - u))), 



-^2^Ai 

{rs^ - in + r2W + (2y - l)?-3 + 2(ri + r2)x), 



2d2imw'^^'^NQ , 3 / , N 2 



L2^/ri 

^"""2^3^/2 (2c?2i((ri' - rs' + l)r3 - 2rix)ri3 + d2r:i{ri^ + (7^2 - r-i)ri^ + (2m - 2r27-3)ri + 2r3(y 

^2''l ''3 ' 

2i2ri2V^ 

/2 

,(2d2i(2(7-i + 7-2) - r3)((ri^ - r2" + l)r3 - 2rix)ri 



(2(i2i(r37'i^ - 2xri - r2^r3 + r3)ri^ + dir:i{r-i^ + 2u) + d2?'3(''3^ - rir-i + 2y - 1)), 

2^2^-1^33/2* 

+?-3^(rf2(r3^ - rir3 - 2r2(ri + r2) + 2y - 1) - di(r32 + 2u))), 
2mw'^^^No, , 3 , J ^ 3 ^ , \ 2 



2mw'^'^NQ 



L2ri^y/ri 
2d2imw'^^^No 



(diira'^ + d22(r-3'' - (ri + r2)r3^ + {2y - l)r3 + 2(ri + r2)x)), 
(dura - 2dii{x + y)r3 + d22(-3r2 + 2r3 + 2(r2 - r3){x + y) - ri{2z + 1))), 
(ri(l-2z) + 2r3z), 

(-3r2 + 2r3 + 2(r2 - r3)(x + y) - ri(l + 2z)), 



■{2d2i{rir3iri + r2 + 2r3) + (2(ri + r2) + r3)u)iiri^ - + l)r3 - 2rix)ri' 



2L2n^r3^/^ 

+r-3^(di(ri(7'i + r2) + u)(r32 + 2u) + d2{r3ri^ + {r2 - r3)r3ri^ 
+ (-37-2r3^ + ur3 - 2r2u)ri - {2r2^ + r3^ + 1)m + 2(r32 + u)y))), 

^^^^^^^f^{2d2ir^{{T^ - r2^ + l)r3 - 2ri.T)(2(ri + r2)x + r3(j/ - z)) - r32(d2(2a;r22 

2i^27'l 7'3 ' 

+a; - 2a;y - r3^{x + 2z) + ri{2r2X - r3{y - z))) - di{r3^ + 2u)(y ~ z))), 

^g^|^(2d2i((ri2 - 7-2' + l)r3 - 2nx){{r3 - 2r2)y + 2r,{x - z))r^ + r3^{d,{r3^ + 2u)y 

-d2i-2y^ + {rs^ - nrs + 2(ri - r2)r2 - 4:x)y + y + 2u + ri{ri - r2 + 4r22:)))), 

^ ANamw'"^r2 
d22 L2^/ri 

AdiiN^mw^/'^J^ 



L2 

(''I +^2), 



d22 L2^ 

2d2iNo7nw'^' ^^^^3 ^ (27^2 - r3)n' + (^2' - r3r2 + 2{u - z))n ~ 2{r3{u - y) + r2y)), 
L2^/ni 
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^io = " ^"t'^^l-^° (c/iir3(ri(ri + ra) + 2u) + ^21(2^32 + r2{2y - 1) + ri(l - 2z)) 

-d22{ri^ + (2r2 - r^y^'^ + {r^^ - rgrj + 2{u - v))r^ - 2(r3(u - y) + rsy))), (178) 

= ^^"'"7"'''^ (.3(-n + .2 + .3) + 2, - 1), (179) 
-L2 

^3, ^ 2mH.V^jVov/Fi ^^^^^^, ^ d22(r3(-ri + + 7-3) + 2y - 1)), (180) 
^2 



A^3 = — — ^— (diir3(2a; + 4y - 1) + d22(-2xri + ri + ra - 2r3 + 2r3X - 2(ri + r2 - r3)?;)), (181) 



^,^ 4iVodnm^V.^ ^^_^^^ ^^^^^ 

I 

Non-zero coefficients and for |(cQ')i['^P,/]) arc 

I 



Al = - ^^^^'^^' ((ri - r2)(r3(2rf2i(ri2 - r22 + l)ri2 + d2(r32 - ir.r^^ + 2ri(ri + r2))) 



+(-4(i2iri3 + 2(i2ri - da^-a)^;) + di(8r3ri3 + 2(r32 - 4r2r3 + 2^)^^ 

-(Srg^ + 2r2r32 + 2wr3 - 4a;r3 + Ar2u)ri + r2r3(r32 + 2u))), (183) 
^2 = ^ — TT^{ri{dir3{rz^ - 2rir3^ - 12ri^r3 + 4rir2r3 - 6Mr3 + 4yr3 - 4riu) 

^1 '3 

-2d2iri((ri2 - r2' + l)r3 - 2rix)(ri(2ri - 2r2 + r3) + 2{u - y))) + d2(4r3ri4 
-2(7-32 + 4r2r3 -2u + 2x)ri^ - {ir^i^ - \r2rz^ + 2(-2r22 + u - 3a; + y + l)r3 

+4r2u)ri2 + r3(r33 + (-2r22 - 2u - 2a; + 3y + l)r3 + 2r2(u + y))ri - r2r32y)), (184) 

or 5/2 

A\ = — (diir3(r2U + riw) + d2i(ri - r2){{rx^ - r2^ + l)r3 - 2rix) + d22(-2r3ri3 



+ {r^^ + 4r2r3 - u - x)ri^ + (r33 - r2r32 + (-2r22 - u - 2x + y + l)r^ + r2(2u + x))rl 
-r2^u + r2r3(w - y) + 2(r32 - .x)(u - y))), (185) 

^2 ^ Sd2iL2mw / („3^^^^3 ^ (^^2 _^ 4^^^^ _ ^ _^ ^y^2 _^ (_r3^22 + (_,,32 + 2m + x)r2 



+r3(r32 - u - 22; + y))ri - r2'^u + r2r^{u - y) + 2(r32 - x){u - y)), (186) 
Aj = -?^5j^^(d2(ri - r2)(2ri - r^) + Admrs), (187) 



^2 ^ ^ ^^^'^^'^\ dnrirs + rf22(ri(r2 - 3ri) - 4u + 4y)), (188) 



Sd2iL2_rnw^^'^ . 
7^ 



Al^ (Sri^ -r2ri +4u-4y), (189) 



A\, ^ ^^^ih^^ir.r, - [r2 + 2r,)r, + 2x), (190) 



^3 

,9 4d2i2"lW^^2 

Al, = V ' (^9^) 
2 2d2-?^2mM/5/2^ 

^12 = 2 ^ — (''I -'"2), (192) 

Al, = ^^^"'^'f (4d2i((ri2 - r22 + l)r3 - 2nx)n' + ^2(^1 - 3r2)r32), (193) 

^1 '3 



-{rirs - (r2 + 2r3)r3 + 2x), 



(rfii(r2 - rijrs + d22(?'i - ?'3?'i - 7'2 + r-2r3 + 2.T - 2)), 



(d2?-3(?-i^ - (r-2 + r3)ri + 2r2r3) - 2d2W(ir^^ - r2^ + l)r3 - 2rix)), 



7* -j^ 2^^3/2 

.5/2 

(ri -r3ri-r2 + r2r3 + 2a- - 2), 



8(i2li2mvi/^/^ , 2 2 



Vr3 



(2d2i(rir3 + 2u)((ri - r2 + l)r3 - 2rix)ri + rs (rir3((i2ri + dir^,) 



+(2diri+d2(ri-3r2))w)), 
2L TO 

' (2d2i(rir3 - (r2 + 2r3)r3 + 2x){{ri^ - + l)r3 - 2rix)ri^ 

^1 '3 

+r3^((i2r3(ri^ - (r2 + r3)ri + 3r2r3) + di{ri - r2)(r3^ + 2u) + d2{ri - 3r2)x)), 
V2' (2d2i((ri" - r22 + l)r3 - 2rix)(ri(2ri - 2r2 + r,) + Au - 2y)r^ 

^1 '3 

+r3^(diri(r32 + 2u) + d2((r-3 - 2r2)ri'^ + (2r2^ - r32 + a; - z)ri + r2(y ~ 4it)))), 

^22 ^21 Cf22 \Ai 

(fi^ - r2ri - 2w), 



8d2iL2mw^^^ , 2 



— {d2i{ri - r2)(ri + r2 - r^i) + durir^ + d22{ri{r2 - 3ri) - 2u + 2y)), 

4^ = J^=8L2m^5/2V?;(ri-r2), 

"22 ^"11 

8L2mw .{^^^^(^j.^ _ r2)(ri + r2 - ?'3) + dii{r2 - 3ri)r3), 

(2d2i((ri' - ra" + l)r3 - 2rix)(y - z)ri^ + di(4r3(3x + 4y - 2)ri'' 



2Liri^^' 

+2((y - z)r32 - 2r2xr3 + 2u{y - z))ri - r3(r32 + 2u){y - z)) + d2(2r-3(y - z)ri'^ 
-((3a: + 6y - ^y^^'^ + 2r2(y - z)r3 + 2x- Axy)ri + r3((x + 2y - 2)r32 + x - 2a;?;))), 

iVomvy"/^ 2 _ ^^2 ^ _ 2ria;)((2r2 + r3)y - 2riz)ri2 + dir3(4r3(2a; + 3y - l)ri2 

+2(r32 + 2r2r3 + 2u)yri - r3(r32 + 2u)y) + d2(4a:ri'' - 2(2r2X + r3(5x + 2y - 'i))ri^ 

+ {{&x + 2y~ 5)r32 + 2r2(a; - 4y + l)r3 + 4u - 8u2;)ri2 + r3(-4?;r2^ + r3(4?/ - l)r2 + (3r32 - 4y + 2)?; 

+&u{2y - l))ri - r32((-2r22 + r^^ - 2y + l)y + u{iy - 2)))), 

— — — (d2i((ri^ - 7'2^ + l)r3 - 2ria;)(?; - z) + diirz{~yrz^ + u + 7-i(ri - r2)a; - 2uy) 

+d22{xri^ + ^3(2/ + 3z)ri^ + {~xr2^ - r3(32/ - z)r2 + u-2{x- l)x - 2(m + x)y - r^'^iy + 2z))ri 
-2^-3^^ + r2U + r-iy - ((ra - r-3)(r3^ + 2r2r3 + 2u) + 2{r2 + r3)x)y)), 

2d2\TllW^ / A / o\3/ri /2 r»\\2/3 2 

z 1= [rz^i - (''2?'3 + 3a;)ri + (2r2a; + r3(-r2 + u + x - 2z))ri + [rzr2 + xr2 

+r3{x + 4y - 2)r2 + 2a;^ - u - 2ux - 2x + 2{u + x)y + r^'^iy + 2z))ri + 2r^y'^ - (rj + (ra^ - l)r3)u 
-(?'3^ + ?'2?'3^ + (-2r2^ + 2u - 2a; + l)r3 - 2r2(w + x))y), 

_ iVoTOw^^/^ 2 - ra^ + l)r3 - 2ria;)ri2 + di{Ar^r^'' + 2(r32 + 2r2r3 + 2u)ri - r^irs^ + 2u)) 

2Lxrx^y/ri 



+d2(2r3ri2 + (-Sra^ + 2r2rz - 4y + 2)ri + rglrs' + 2y - 1))), 

2d2imw'^/'^No 2 2 , IN o ^ 
z — ^ {{ri ~r2 + 1)7-3 -2rix), 

-diir3{ri^ - r2ri + 2m)) 
2d2imw'^^'^No , 3 , ,2 



(d22(?'i^ - (2r2 + r^yi^ + (r2(r2 + rg) + 2{u + x - l))ri - 2(r2 + rg)^ + 2(2r2 + r3)y) 
•i+2m)), 

(ri^ - (2r2 + r3)ri^ + (r2(r2 + rg) + 2{u + x - l))ri - 2(r2 + r3)u + 2(2r2 + ra)^), 



(?-3^ + (?-2 - ri)r3^ + (2y - l)r3 + 2(ri - r2)a;), 



2d2imw'^^'^No , 3 , ,2 



{2d2i{{ri^ - r2^ + 1)^3 - 2rix)ri^ + da^alri^ - (r2 + r3)ri2 + 2(r2r3 + u)ri + 2r3{y - u) 



{2d2i{r3ri^ ~ 2xri - r2^r3 + r3)ri^ + dirgir^'^ + 2u) + d2r^{r^'^ - rirs + 2y - 1)), 



(2d2i(2ri - 2r2 - r3){{ri' - r2^ + l)r3 - 2rix)ri-' + r3^(d2(r3^ - nrg 



2Liri2r33/2 
+2(ri - r2)r2 + 2y - 1) - diir^^ + 2u))) 
2mw^^^No 



2mw'^''^NQ 

Nod2mw'^''^ 
2d2imw''^^Nn 



[diir^i'' + rf22(r3'' + (r2 - ri)r:i' + {2y - l)r3 + 2(ri - r2)x)), 



(dnr3(l - 2z) + d22(-3r2 - 2r3 + n{l + 2z) + 2(r2 + rz){x + y))), 
(ri(l-2z) + 2r3z), 

(-37-2 - 2r3 + ri(l + 2z) + 2(r2 + n){x + y)), 



■(2d2i(rir-3(ri - + 2r3) + (2ri - 2r2 + r:i)u){{ri'^ - r2^ + l)r3 - 2rix)r{'' 



2Liri2r33/2 

+r3^(di(ri2 - r2ri + u)(r32 + 2u) + d2{rzri^ - r3(r2 + rz)ri- 
+(3r2r32 + ur3 + 2r2u)ri - (2r2^ + + l)u + 2{r-^^ + u)y))), 

N^mw ' J rz'{d2{2xr2' +x- 2xy - r^^^x + 2z) - r^{2r2X + r^{y - z))) 
2L\r\^r-};^i ^ 

-diir^'' + 2u){y - z)) - 2d2in\{n^ - r2^ + l)r3 - 2rix)(2{n - r2)x + n{y - z))), 

^omw ' j ^^2(^^_^^_2y2 _ (_^g2 _^ ^^^^ _^ 2r2(ri + ra) + 4x)2/ + y + 2^^ + ri(ri + ra - 4r2x)) 

-rfi(r32 + 2u)y) - 2d2iri'((ri2 - r2' + l)r3 - 2ria:)((2r2 + r^)y + 2ri(a; - z))), 
_ 4iVor2mty^/2 

^22 -^iV^ 

4diijVomvy^/2^ 

- -('-1 - r2), 



d22 Liy/r^ 

(ri^ - (2r2 + rs,)ri' + (r2(r2 + rg) + 2( u — z))ri — 2r3U + 2(r2 + r^)y), 



2d2iNomw'^^^ 3 , ^ 2 



^ ° (diir3(ri2 - r2ri + 2u) - rf22(ri^ - (2r2 + r3)ri2 + (r2(r2 + r^) + 2{u + x- l))ri - 2r3M) 
-2d22{ri + ra + rj,)y + d2i{-2yr2 + r2 + 2r3Z + ri(l - 2z))), 
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= p — +?'2 - rsjra - 2y + 1), (230) 

^1 

4; = - ^"'^'^'^"^ (dnrs^ + d22(r3(-ri - + rg) + 2y - 1)), (231) 

OlSl m 

= ° )L (dura (-2a; - 4y + 1) + d22(-22/r2 + + ir^z + ri(l - 2z))), (232) 

A|;^- ''^"^°"r'''^ (,-z). (233) 
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